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1 Introduction 



A Ricci flow on a manifold M is a family of Riemannian metrics gij{t;x), x £ M, 
t £ I QM the family parameter, and / a connected interval, satisfying 

^ = -2i?., , (1.1) 

where Rij is the Ricci tensor of gij{t,x). It is usually more convenient to study the 
Hamilton-DeTurck flow (or Ricci- DeTurck flow) 

^ = -2i?,, + £xg^J . (1.2) 

If g solves (jl.2p . then the pullback ip^g solves (jl.ip ([2j, p 80), where ^pt is a family 
of time-dependent diffeomorphisms generated by the vector field X. From (|1.2p . the 
scalar curvature R of gij{t) evolves according to 

— = AR + VxR + 2RijR'^ , (1.3) 

where A := g'^^ViV j is the (t-dependent) Laplacian. 

Ricci flow of asymptotically flat manifolds arises in several physical and mathe- 
matical contexts, ranging from the physics of closed string tachyon condensation jS] 
to existence problems for static Einstein metrics [7J. It is known that asymptotically 
flat initial data remain asymptotically flat and smooth when evolved by Ricci flow for 
some time interval [0,T) [11[I3]. For complete, rotationally symmetric initial data with 
no minimal sphere present (corresponding in static general relativity to the absence 
of black hole horizons), it has been shown that the flow exists for all future time and 
converges to flat space [13] . 

Then the question arises as to what happens when complete, rotationally symmetric 
initial data containing a minimal sphere are evolved. There are two cases, depending 
on whether the minimal sphere is topologically essential (representing a nontrivial class 
in 7r2(M)) or, as will be the case for data on M", inessential. Of these, the essential 
case is the easier one to study, and is the subject of the present paper. 

It is well-known that a closed 3-manifold admitting an essential minimal 2-sphere, 
when evolved under the Ricci flow, will develop a singularity within finite time. Depend- 
ing on the initial configuration, the singularity may be localized at or near the minimal 
sphere, or may be global in the sense that the manifold collapses everywhere at that 
time. It is expected that this result carries over to the noncompact case. Somewhat to 
the contrary though, asymptotically flat initial data are critical data with respect to 
the pinching theorem of M Simon [T3]. That is, Simon proves that warped products of 
a line with a positively curved closed manifold will evolve to form a neckpinch singu- 
larity in finite time, provided the initial data obey certain asymptotic conditions at the 
ends of the line, including a condition on asymptotic growth of the area of the closed 
manifold factor. Simon's growth condition is written as a strict (i.e., open) inequality 
which is not satisfied by asymptotically flat manifolds, but such manifolds would lie in 
the closure. This suggests that such data may exhibit interesting evolutions, includ- 
ing perhaps critical phenomena of the sort observed in certain numerical evolutions in 
general relativity in a scenario sometimes called "critical collapse" jl]. 
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Garfinkle and Isenberg [S] studied numerical Hamilton-DeTurck flow of a 3-sphere 
with "corsetted" initial metric, admitting an inessential minimal 2-sphere (the "waist"). 
For tight corsetting, meaning that the waist has very small area relative to the 2/3'''^'^- 
power of the volume of the 3-sphere, the waist "pinches off" (a local singularity forms 
there), whereas for more gentle corsetting, the entire sphere shrinks to a point before 
a local singularity can form. The critical solution separating these two alternatives is 
a degenerate neckpinch singularity modelled by the Bryant soliton [6]. 

The Garfinkle-Isenberg result shows that initial data for the Ricci flow with an 
inessential minimal sphere divides into two disjoint sets whose common boundary con- 
sists of points whose evolutions exhibit critical behaviour. 

Husain and Seahra [10] then considered a numerical Hamilton-DeTurck flow of a se- 
quence of initial metrics, each rotationally symmetric and reflection symmetric through 
an essential minimal 2-sphere, sometimes referred to as the "throat" or "bridge" . Evo- 
lution occurred on a bounded region, with boundary conditions imposed. Two alter- 
natives were again found; for some initial data, the throat pinched off while, for other 
data, the throat expanded to infinity. 

We study a related question. Consider simply connected, noncompact, complete 
manifolds with rotational symmetry. These are M", M x S"'~^, and certain quotients 
thereof. As M" has no essential minimal sphere, consider M x and endow it with 

an SO(n)-symmetric metric such that there is a minimal hypersphere located at, say, 
r = and an isometry corresponding to reflection in that hypersphere. One such metric 
is the t = slice of the (n -|- l)-dimensional Schwarzschild-Tangherlini metric. One can 
then identify points under the action M x S^~^ B {r,p) i— )• (— r, —p), where p G S^~^ and 
—p is the antipode of p. This produces a smooth metric on M x MP""-*^ = MP" \ {pt}. 
For n = 3, this model is known in gravitational physics as the MP^ geon. We pose and 
will answer the question, "What is the Ricci flow evolution of the MP^ geon?"lf| 

In section 2, we adapt a standard Ricci flow argument from the setting of closed 
manifolds to that of asymptotically flat manifolds. This forms the basis for what 
follows, and shows that essential 2-spheres in asymptotically flat 3- manifolds collapse 
whenever the initial scalar curvature is bounded below by a nonpositive constant that 
can depend on the initial area of the minimal sphere. The result is valid for arbitrary 
Hamilton-DeTurck flow, including Ricci flow. As the covering space of the MP^ geon 
obeys this bound, this suffices to answer the question just posed. 

This raises the possibility that, by choosing initial data for which the lower bound 
on initial scalar curvature is violated (e.g., by choosing a different initial metric on the 
geon manifold, with the same isometries), interesting dynamics might arise, such as 
observed in the studies cited above [U EJ [10] . It also raises the question of the precise 
comparison, if any, of this result to the numerical work of [lOJ, who found that for 
certain initial data with scalar curvature well above our lower bound, their numerical 
Hamilton-DeTurck evolution did not lead to collapse. 

To understand these issues, we perform our own numerical simulations. In section 
3, we lay the groundwork. We first discuss Hamilton-DeTurck flow with the "DeTurck 
trick" formulation which we use for our numerical evolutions in Section 4 (and which 
[5] used). For comparison purposes, we then discuss the normal coordinate Hamilton- 

^This was posed by DM Witt to the second author quite some time ago. 
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DeTurck formulation used in |10j . 

Section 4 contains our numerical results. We use the same form of initial data as 
Husain and Seahra [lOj and, for practical purposes, we also now restrict our evolution to 
a bounded manifold, but we use different evolution equation^ and different boundary 
conditions at the boundary at large r. Our numerical evolutions always exhibit collapse 
of the throat. Section 5 contains a brief discussion of the numerical evolution. 

Our convention for the curvature tensor R"'hcd is that used in |9j and equals the 
quantity denoted by Rcdb^ in [2j. We write Rabcd '■= daeR^hcd- We denote the Laplacian 
by A := g°'^W gS/h- For a definition of asymptotic flatness, see |13j . 



2 Time derivative of the area of a minimal sphere 
2.1 The maximum principle for scalar curvature 

For Ricci flow on compact manifolds, it is standard that positive scalar curvature is 
associated to contraction and concentration of curvature under the Ricci flow. In this 
subsection is that scalar curvature, we recall that this is also true for complete manifolds 
with vanishing scalar curvature at infinity. 

Proposition 2.1. Assume that (II. 3|) has a solution on some time interval [0,r], 
T > 0, on a complete manifold M with one or more asymptotic ends, and assume that 
R{t,x) — >■ as X ^ oo (that is, as the point x tends to an asymptotic end), for all 
t G [0, r]. Let the scalar curvature Ro{x) of the initial metric goix) := g{0,x) obey 



inf {Ro{x)} =: -a^ , (2.1) 



for some a > 0. Then 



-a" 



^(*'^)^^^^- (2.2) 
1 + —t 



Remark 2.2. The existence assumption always holds for Ricci flow developing from 
asymptotically flat initial data ([13]. [3]). 



Proof. Equation (jl.Sp can be written as 

1 2 

+ VxR- (2.3) 



dR ^ „ 2 1 2 



AR+-R^ + 2 
at n 



Rij ^ gij R 



Exhaust M by a sequence of closed, bounded sets Ki. An easy application of the 
maximum principle (e.g., [11]) to the i^^ subset yields that the minimum of R occurs 
on the parabolic boundary of [0, T] x K^. Now take i — )• oo and use that R{t, x) — ?• as 
X — 7> oo to deduce that the minimum must either be zero or must occur on the initial 
data. This proves that 

R{t,x) > inf {Ro{x)} = -a^ a > 0. (2.4) 



^We use a difTerent Hamilton-DeTurck system. Occurrence of collapse of the minimal surface should 
not depend on this choice; see section 2. 
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If a = 0, this proves the theorem, so we now consider a > 0, in which case the infimum 
is a minimum. 

For a > 0, let Q := — ^ ^1 + R- Since R = —a^ < at its minimum, then 
Q > there. Hence the maximum of Q is positive, and it then follows from the 
definition of Q that, at the point where Q achieves its maximum, R is negative. As 
well, Q{t, x) ^ as X — )• oo. 

Work on the intervalcompact set [0,r] x Ki. From (j2.3p . we have 

^<^Q + VxQ + -R{Q-l) ■ (2.5) 
ot n 

Then from the maximum principle (and since R < where Q achieves its maximum), 
either the maximum of Q occurs on the parabolic boundary of [0,T]x Ki or Q — 1 < at 
the maximum. If it occurs on the parabolic boundary, by taking i ^ oo and using that 
Q{t,x) — )• as x ^ oo, we see that the maximum must occur on the initial boundary, 
where Q{0,x) = —-^R{0,x) < 1. Hence, Q{t,x) < 1, and the result follows. □ 

Remark 2.3. Obviously, these results also hold on [0,T] x K , K C M a compact 
set with boundary dK, provided that R(t, x) = for x G dK. But they do not hold 
if merely R{t, x) > —const on dK, an observation which appears to be relevant in 
numerical studies (see section 3.2.2). 



2.2 Evolution of hypersurfaces 

In this subsection, we obtain the evolution equation for the area of a closed hypersurface 
T, in an n-manifold {A4'^,g), defined (locally at least) by an expression of the form 
F(x*) = 0, where are the coordinates held constant in the derivative ^ appearing 
in (fr2D|5| 

Let be a smooth unit vector field normal to S. Let hat be the induced metric and 
let Hab be the extrinsic curvature of S, with trace H (thus H is the mean curvature, 
taken here to mean the sum rather than the average of the principal curvatures at a 
point). If Hab denotes the trace-free part of Hab, then 

Hab = Hab H -rhabH . (2.6) 

n — 1 

If TZ denotes the scalar curvature of the induced metric hab on S, then 
R — Ric(n, n) = — 

by the Gauss equation. Finally, let the induced area element on S be dfi. From the 

''For hypersurfaces defined more generally by an equation of the form F{t, a;*) = 0, the generalization 
of our result is easy to obtain, using the transport theorem. 



n 



n 



H^ + R+ \H\ 



(2.7) 
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Hamilton-DeTurck flow equation (jl.2p . the area |E| := j-^dji of S evolves as 



' ' hdt ^ 2 h dt ^ 2 L dt ^ 



- [ h'^ {Rij - ViXj) dfi = - [ {R- Rijn'n^ - h^^ViXj) dfi 



using ()2.7p . the divergence theorem on S, and -ff := h^^Vinj. 
2.3 3-manifolds 

Now fix the dimension to be n = 3. Then we have J^TZd^ = 47rx(S), where x(^) is 
the Euler characteristic of S. Then 

= -2^X(S) + - ^ / I^p) / Hn-Xd^i, (2.9) 

where W{T,) := \ H^dfj, is the Willmore energy [151 of S@ Using Proposition 2.1, 
()2.9p becomes 

<-2^x(S) + i^(S) + -^^ + J^Hn-Xdfi. (2.10) 

The tightest bound occurs in the case of a minimal 2-sphere. Then H\y and the 
Willmore energy vanish and xi^'^) = 2 so 

d , , 3a2 ^ ^ 

— S <-47rH S . 2.11) 

We note that the initial time derivative of |S| will be negative whenever 

n Stt 

a' < ^ , (2.12) 
l^lo 

where \Ti\q denotes the initial area of S. Even when not initially negative, the derivative 
can turn negative and remain so thereafter. 

Now we can regard the left-hand side of (j2.1ip as the forward difference quotient of 
the functional whose value is the area of the smallest essential minimal surface present 
at t. Then (see Lemma 2.22 of |12j ) an upper barrier for the area of the smallest 
minimal surface at any t > is provided by the function ^(t) that solves the initial 
value problem 

d'^ Sa^^- 

= -47r , 2.13 

dt 6 + At ' ^ ^ 

^(0) = |S|o . 



When E is a 2-sphere, the combination of Willmore energy and x(^) that appears in (|2.9p also 

|S|l/^ 
167r3/2 



appears in the definition of the Hawking quasi- local mass niH ' . [47r — Vy(S)] 
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The solution is 

^(,) = I #^(1 + 1^) + (l^lo - sfe) (1 + ¥r^" for / 4/3, 
\(1 + ft) [|S|o - evrlog (1 + ft)] for = 4/3. 

Proposition 2.4. Say that the conditions of Proposition 2.1 hold and that {M^go) 
contains a minimal sphere of area 4tt6 for some 6 > 0. Let 3a^ — 4 < 6/6. Then the 
Eow fails to exist in finite time. 

Proof. We must show that if the flow lasts long enough, the minimal sphere collapses. 
If a-^ = 4/3, this is immediate from the second line of ()2.14p . since ^ has a zero at some 
t > 0. Otherwise, let 7 := Sa^ - 4 / 0. The top hue of (pT4]l becomes 



^{t) = 47r ( 1 + 



so ^(t) = for some t > iff 



2 y/' 6 



(2.15) 



Now >lif0<7<6/(5 and < < 1 if 7 < 0. The left-hand side equals 1 at 
t = and is otherwise monotonic, increasing without bound if 7 > 0, and decreasing 
with asymptote if 7 < 0. Hence there is always a root for t > if 7 < 6/5. □ 

Corollary 2.5. For any X, consider the flow (jl.3p developing from an initial metric 
which describes the RP^-geon. Then the flow fails to exist in finite time. 

Proof. The Riemannian double cover of the geon is the 5 = 1, a = case above. 
Furthermore, since the initial data are asymptotically flat, then R{t,r) as r ^ 00 
for all t during the flow ([HUB]). Then the result holds on the Riemannian double cover 
and, hence, on the manifold itself because the covering map is a local isometry. 

Alternatively, the geon itself has xi^) = 1 and |S| = 27r, so equation (2.13) 
applies, with the 47r on the right-hand side of the differential equation replaced by 27r 
and |So| = 27r. The solution (|2.14p is then multiplied by an overall factor of 1/2, which 
has no effect in the proposition. □ 



2.4 2-dimensions 

The n = 2 case is itself of some interest. Then (12.81) becomes 



dt 



S| = -- /" {R-2Hn- X)dn 



(2.17) 



with H the geodesic curvature of the closed curve S. Then 

dt 2 JY, 



'S| = — / Rd^L (2.18) 
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if S is a geodesic. If i? ^ at infinity (or on any boundary), Proposition 2.1 yields 

2 

that R > —jf^, and so (2.19) becomes 

d , , , , , , 

— S < ^ S . 2.19 

dr ' - 2(1 + a2t)' ' ^ ^ 

Thus 



\T.\{t) < m{0)Vl + aH . (2.20) 



3 Rotational symmetry and Hamilton-DeTurck flow 

There are two scales in this problem, the initial area of the minimal surface and the 
minimum of the initial scalar curvatureQ The metric can be rescaled to fix one but not 
both of these. By choosing 5 = 1, we fix the initial area. There is then the intriguing 
possibility that collapse of the minimal sphere will not occur if the geon initial data is 
replaced with initial data with scalar curvature sufficiently negative that > 10/3. In 
the next section we will study this question numerically. In this section, we formulate 
the equations we will need. 

There are several versions of Hamilton-DeTurck system that are used in the study 
of rotationally symmetric Ricci flow. We discuss the formulation we will use for our 
numerical integrations in the first subsection. We then discuss, for comparison purposes, 
other formulations that have been used in related work. 

3.1 DeTurck's background connection method 

3.1.1 The flow equations with a background connection 

The first Hamilton-DeTurck system arises from DeTurck's original trick for proving 
short-time existence of Ricci flow, and results in a parabolic strictly system of two 
equations. The idea is that parabolicity fails only because there are families of flows 
whose members are distinguished from each other only by a continuous time-dependent 
deformation of the coordinates. Thus, parabolicity is restored by "breaking coordinate 
invariance" . This is done by fixing a t-independent background connection throughout 
the flow. That is, let F*^ be the flowing connection (compatible with gij{t)) written 
in local coordinates, and let f*-^ denote the chosen background connection in these 
coordinates. Define the vector field 

X' := g^' (rl^ - f),) . (3.1) 

Then the system (|1.2p is parabolic (for any choice of t-independent connection T). For 
rotationally symmetric flow on M", an obvious choice of coordinates is 

ds^ = e2^(*'^)dr2 + r2e2^(*'^)5(5"-\ can) , (3.2) 

where g{S"'~^ , can) denotes the metric of unit sectional curvature on S"""^. Then it 
makes sense to choose a background connection to arise from a rotationally symmetric 

^This presumes the scalar curvature has a negative minimum. In the initial data we study, the only 
exception will be the Schwarzschild-Tangherlini data, for which R(0,p) = for all points p. 
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background metric. Since every rotationally symmetric metric is conformal to a flat 
metric, we write the background metric as 

g = e2^('-) [dr^ + r25(5"-i, can)] , (3.3) 

for some function il){r) that must be specified. Choices of ip with rotational symmetry 
and asymptotic flatness include 



, 



flat background, 



|ln[(l +r"")] , a = const>Q. 



(3.4) 



The last choice listed is a family of metrics which were used as initial data for Ricci 
flow (with 1 < a < 2 and n = 3 dimensions) in [10] and which would, if a = n — 2, 
correspond to a time-symmetric slice of (n + l)-dimensional Schwarzschild spacetime. 
Then the DeTurck vector fleld is 



^ ' ^ dr 



-2A 



dA 

57 - ^ (0 



(3.5) 



n 



g-2A _ g-2B 



+ e-2^^-e-2V(r) 
or 



d_ 

dr 



The Hamilton-DeTurck system becomes 



dA 
'dt 



dB 
'dt 



-2A 



d^A 

Qj.2 



dA 
dr 



+ (n-l) 



dB 1 

dr r 



, OR ri9^ 2dB 1 ,,,, , dA , dB ^ 
+ n - 1 2 + V-^- + - 2— V^' r 



-2A 



d^B 

Qj.2 



+ (n - l)e 
If we deflne 



-2B 



/dB 1 

\ dr r 



1 



{n-2) ^_2B 



(3.6) 

(3.7) 
(3.8) 



to be the area of r = const spheres, with Vn-i being the volume of an (n — l)-sphere 
with unit sectional curvature, we can rewrite the above system with B replaced by S. 
In n = 3 dimensions, this gives 

'dAy ^ 1 fdSV 
dr 



dA 
'dt 



-2A 



d^A 



+ - 



S 



Qj.2 

IdA 

r dr 



^ 252 ^dr) 



dA 

V"(r) + ^V'(r) 



1 dS 1 , 
rb dr 



dA 2 
dr r 



i^'{r) 



1 dS 



V''(r) 



dS_ 
dt 



-2A 



S dr 

\ d^s 

I Qj.2 S \dr J 



I /dS\' 



dS I 87rr2 dS 



5 dr 



(^\r) + - 
\ r 



(3.9) 
871(3.10) 



Once a function ^ is chosen, it is clear that this system is parabolic. This system was 
used in [5] in their study of rotationally symmetric flow on "corsetted" 3-spheres. 
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3.1.2 Minimal surface boundary conditions with background connection 



Let us now consider the problem of the previous section, which is the evolution of an 
asymptotically flat metric on a manifold with inner boundary T, which is a minimal 
surface. At S, both the mean curvature and the vector field X must vanish, so 



:= Vit,ro) 
H^{t) := H{t,r^) 



The second of these conditions is of course ^ 



= , 
= . 

or, equivalently. 



dB 
dr 



+ 



ro 



, 



while the first then yields 



-2A{t,ro) 



dA 
dr 



i^'iro) ) + (n - l)e-2^(*''^'') ( ^ + V'(ro) 1 = 



(3.11) 
(3.12) 



(3.13) 



(3.14) 



3.2 Polar coordinate gauge 

3.2.1 The Row equations in Gaussian polar coordinates 

This system was used by Husain and Seahra [TO] and is not parabolic. The metric is 

ds^ = dr^ + F^{t,r)g{S"-\ can) , (3.15) 

corresponding to A = 0, B = log{F/r) in (j3.2p . The coordinate r is distance from 
the minimal sphere at r = 0. In this system, the Hamilton-DeTurck fiow becomes a 
constrained system comprised of the differential equation 



dF d^F 1 fdF 



+ — 



dt dr"^ F \ dr 
and a differential constraint 



V- 



dF n-2 



dr 



F 



dV 

dr 



(n - 1) d^'F 
F dr'^ 



(3.16) 



(3.17) 



which determines the generator X = V-^ of the diffeomorphism in ()1.2p . The mean 
curvature of a constant-r sphere is 



H :-- 



{n-l)dF _ 1 dS 
F 9^ ~ S~d^ 



(3.18) 



where S = Vn-iF^ ^ is the area of an r = const sphere (cf (j3.8p ). and so the solution 
of the constraint can be written as 



r 

V{t, r) = C{t) + H{t, ro) - H{t, r) - [ H\t, r')dr' , 

[n-l) J 



(3.19) 
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where C{t) is an arbitrary function of t, to be determined by boundary conditions at 
tq. When the manifold contains an origin for the rotational symmetry (e.g., M"), this 
formulation suffers from the problem that H diverges there. 

In the case of n = 3 dimensions, (|3.16p leads to the remarkably simple equation 

dS__^ ^dS_ _ 
dt dr"^ dr 

for the area S{t, r) of an orbit of the rotational symmetry. Indeed, this equation would 
be linear if not for the dependence of y on 5* via the constraint (j3.17p (which takes a 
slightly more complicated form when F is replaced by S). 



3.2.2 Minimal surface boundary conditions in polar coordinates 

The condition that V should vanish at the boundary S is simply C{t) = 0. The 
condition (j3.12p that the boundary S be a minimal surface is H{t,rQ) = 0. These 
conditions, applied to (13.19p . then imply that 

r 

V{t, r) = -H{t, r) - I H\t, r')dr' . (3.21) 

ro 

Since H is not identically zero (for all r), then V[t,r) < on every convex surface 
r > tq. Husain-Seahra |10] impose that F — t- at infinity and then approximate this 
for numerical purposes by imposing V{t, Tc) = on an outer boundary r = Tc at all 
As a result, they have 

Tc 

Hit, r,) = -^-^ [ H\t, r)dr < , (3.22) 
[n-l) J 

and = only when H{t, r) = for all r, so the r = Vc boundary is necessarily concave 
in the direction of increasing r at all t > 0. This has the potential to generate negative 
scalar curvature at r = rc since 

R = -2^ - r^^H' + (n - l)(n - 2) (S/K-i)-'/^""') • (3.23) 
or [n — 1) 

Then Proposition 2.1 and, concomitantly, Proposition 2.4 would no longer apply. 



3.3 Area radius coordinates 

This is also not a parabolic system, but can be reduced to a single parabolic equation. 
The technicque was used by [T3] to prove existence and convergence of rotationally 
symmetric asymptoticall flat flows. Unfortunately the technique fails for initial data 
containing a minimal surface. For this reason, the results of |13j do not apply when 
minimal surfaces are present. In the rotationally symmetric setting, fix X such that 
the components of the flow equation in directions tangent to the orbits of rotational 

^In [To] the outer boundary is defined by F{t, rc) = const and thus rc can vary in time. 
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symmetry are trivial; i.e., so that = for i, j 7^ 1. This means that the metric 
coefficient B above wiU be time-independent. We let f = and so we write the metric 
as 

ds'^ = f {t,r)dr^ +r^g{S''-\ can) , (3.24) 

and call r the area radius. Then in local coordinates this fixes 

R22 d _ d _ 

- guTl, dr - guFl, dr ' ' ' ' ' ^'-"'^ 

Then (|1.2p becomes a constrained system in which the constraint can be solved, leading 
to the single parabolic equation 

dt /2 9r2 P\drJ ^\ r rp ) dr r^f >' ^ ' 

Having a single parabolic equation instead of a system is an enormous advantage |13j . 
and the maximum principle can be used to show that / remains bounded whenever it 
is bounded on the initial data, yielding uniform parabolicity. However, unboundedness 
of / on the initial data corresponds precisely to the presence of a minimal sphere, since 
the mean curvature of r = const sphere is given in this system hy H = ~p ■ (It is 
shown in [13] that no solution of (j3.26p can form a minimal sphere during the evolution 
if none is present initially.) 



3.4 Conformal gauge: n = 2 dimensions 

Every rotationally symmetric metric is conformally flat. However, for n > 2, the 
conformal class varies throughout the flow. To see this, set A = B in the metric (13.2jl . 
Then equations (j3.6l 13. 7p reduce to the single equation 



dA _ 
dt \ dr"^ 



-2A 



r d'^A {n - 1) dA 



+ 



dr 



+ (n - 2)V''(r) 



dA 1 
dr r 



(3.27) 



together with a restriction on ^: 



= (n - 2) 



^4) +r{r)-24^\r)-U'{r) 



(3.28) 



When n > 2, we obtain from (|3.28p that 



C{t) 



^— (t,r')j dr' 



(3.29) 



for some function C{t). Since the left-hand side is time-independent, so must be the 
right. Thus, the system (j3.27l I3.28P is rarely solvable unless n = 2. In n = 2, however, 
(j3.28p is trivial and then (j3.27p can be solved. 
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4 Numerical results 



4.1 Initial data 

We begin with the class of metrics 

= r^^^ + p'9{s--\^^^) , (4.1) 
1 - [po/pr 

where g{S"'~^ , can) is the metric with constant unit sectional curvature on the (n — 1)- 
sphere. For n = 3 and a £ [2, 1], these metrics were used as initial data by [10]. When 
a = n — 2, this is a static hypersurface in the Schwarzschild-Tangherlini metric. There 
is a minimal surface at p = pQ. 

To obtain initial data for the functions A{t^r) and S{t,r) used in the evolution 
equations, we first transform to isotropic coordinates. We obtain 

ds2 = p\r){dr^ + r^g{S''-\ can)) , (4.2) 

where ro is a constant of integration arising in the coordinate transformation. Areas of 
constant-r spheres are given by 

(^)*(..|,)^, (4.4, 

where Vn-i is the volume of an (n — l)-sphere of constant unit sectional curvature. The 
minimal sphere now lies at 

r := ro/2i/° . (4.5) 

We choose ro = 2^/", so the minimal sphere occurs always at r = 1. We then fix po = 1; 
which fixes the area of the initial minimal sphere to be Vn-i- Then the initial data are 

AiO,r) = log/3 = -logfl + 4') --log2 , (4.6) 



1 / 1 \ (n-l)a / I \ {n-l)a 

S{0,r) = V{r)r--'{-] + ^) ■ (4-7) 



4.2 The boundary data 

Our initial data have reflection symmetry in the minimal sphere, as well as rotational 
symmetry. The flow equations preserve isometrics, so the subsequent evolution will 
share these symmetries. 

We therefore choose to place a boundary at the location r = 1 of the minimal sphere. 
The idea is then to take the Riemannian double {D^gjy) (the topological double, with 
the Riemannian metric qd induced by pullback of g{t) under the covering map) at any 
time t along the evolution. For this, we must keep the mean curvature zero at r = 1. 
Rotational symmetry then ensures the r = 1 sphere is totally geodesic, so the double 
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is a smooth manifold. In general, we do not expect R{t, 1) = at the r = 1 boundary. 
Nonetheless, the arguments of Section 2 apply to {D,g£){t)) provided that R{t,r) — )• 
on the asymptotic ends (D, of course, has two of them). 

As we wish to study the noncompact case, we would like to allow r to range through 
all values > 1. However, for numerical purposes, we must either choose a finite cut-off 
or use a more sophisticated method (such as conformal compactification, which then 
would introduce the complication of dealing with a singular boundary value problem). 
For simplicity, we choose a finite cut-off, so r G [l,?"c]- 

We now have two boundaries and a parabolic system of two PDEs, so we expect 
four boundary conditions, two at each boundary. Three are obvious. These are that 
the mean curvature should vanish at the r = 1 boundary and the DeTurck vector 
field should vanish at both boundaries (to make it possible to infer the appropriate 
conclusions for Ricci flow from our results for the Hamilton-DeTurck flow on the same 
bounded manifold). For the remaining condition at the outer boundary, we would 
prefer to hold the scalar curvature constant (preferably zero) there. Instead, we must 
set a condition that contains no worse than first derivatives of A and S. We choose 
to hold the mean curvature equal to a constant A on the outer boundary. For the 
case of a single boundary, Cortissoz [3] showed long-time existence of the flow with 
this boundary condition, so it would be interesting to see if any vestige of this result 
remains for the present case (speciflcally, noncollapse of the minimal surface when a is 
large) . 

The conditions become, for all t > 0, 



Atr = l: {9r (4.8) 




const 

We take A to have the value that it would have in the background metric at r = Tc, so 

(4.10) 



A = (n - l)e-'^("=) (^^ + i^'irc)^ 



4.3 The numerical results 

For the initial data of section 4.1 and n = 3 dimensions, the initial minimal surface is 
always of area 47r, so 5 = 1 and thus = 10/3 is the critical value in Proposition 2.4. 
This value is achieved for a = 8/3. 

Our flrst numerical run probes the case of n = 3 dimensions, with a taking values 
below or near the critical value. Figure 1 shows the results for a values 1, 2, 2.5, and 
3. In every case, the minimal sphere collapses to zero area. Though we do not show 
it, we monitor the scalar curvature at the position of the minimal sphere and observe 
that it diverges to +oo at the collapse time. 

In the second run, we probe larger values of a, in the hope of seeing expansion of the 
minimal surface. The results are displayed in flgure 2. We use a values all the way up 
to a = 8. In each case, collapse eventually commenced, after which there is no evidence 
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alpha = 1 
- - alpha = 2 

alpha = 2.5 

alpha = 3 




0.2 0.4 0.6 0.8 
Time 



1.2 1.4 



Figure 1: Area of the minimal surface at r = 1 as a function of flow time, for a = 1, 2, 
2.5, and 3. 




Figure 2: Area of the minimal surface at r = 1 as a function of flow time, for a = 4, 5, 
6, and 8. 
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Figure 3: Area of the minimal sm'face in n = 2 dimensions as a function of flow time. 
The initial data are again given by (j4.61 14. 7p with n = 2 and various values of a. 

for subsequent re-expansion. However, at large a, the expansion appears eventually to 
slow and halt. At such values, as the evolution proceeds, sectional curvature becomes 
highly concentrated at the minimal surface, manifested as a large value of the second 
spatial derivative of the area S. This may indicate that a singularity forms before 
collapse of the minimal surface occurs, but our numerics are not sufficiently reliable 
when derivatives become large. 

Our third graph deals with the case of n = 2 dimensions. For initial data derived 
from the form of ()4.ip with a > 0, the initial scalar curvature will always be negative in 
2 dimensions (as is the case in any dimension when a > n — 2, so each curve in figure 3 
will initially expand. By ()2.20p . we do not have a collapsing upper barrier function for 
any a, but it is still possible in principle that collapse could occur. We find no collapse, 
despite running for much longer times than for the higher-dimensional cases. 

The numerical integrations were performed with MATLAB's "pdepe" integration 
routine for partial differential equations. We performed several tests of the validity of 
the code. From equations (jS.lOp . (|4.6p . and (|4.7p . we compute ^(0, 1) = 47r(a — 2). 
This serves to define our unit of time, but also serves as a test of validity of the code, 
in the sense that the initial derivative should be linear in alpha and zero at q = 2. This 
is verified by figure 4, which presents a close-up of the evolution near the initial time. 
For example, at a = 2 the slope is zero initially. 

As additional verification, we vary the dimension, using initial data for which the 
scalar curvature obeys R{0, r) = in each case. The initial data may therefore be 
regarded as the metric on a moment of time symmetry in an (n + l)-dimensional 
Schwarzschild-Tangherlini exterior spacetime. We therefore use the evolution equations 
()3.6H3.8p with arbitrary n > 3 and set a = n — 2 in (14.61 14. 7p . For these runs, we 
normalize the initial area of the minimal surface to be 27r"/^/r(n/2). It is easy to verify, 
using equation ()2.8p . and the fact that R{0,r) = for such data (and thus R{t,r) > 
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Figure 4: Area of the minimal surface at r = 1 as a function of flow time for times 
close to zero, for q = 1, 2, 4, and 8. 



by the maximum principle), that for asymptotically flat boundary conditions the area 
of the minimal surface will be bounded above at all times by 



27r"/2 



--{n-l){n-2)t 



n > 3. 



(4.11) 



This should therefore give an indication of the effect of our boundary condition at finite 
r = Tc, which we typically pick to be at = 10 (where r = 1 is the minimal sphere). 
The results are depicted in figure 5. 

A further technical check concerns the choice of different Hamilton-DeTurck flow. 
For this, we return to (13. 4p and this time choose the flat background, so ip = Q m. 
the evolution equations. The result is displayed in figure 6. The evolutions take much 
longer to collapse, but collapse occurs just as before, and in fact the code seems more 
reliable close to the collapse timejl We present graphs for small a, but we carried out 
the calculation for larger a as well, with similar results. 



5 Discussion 



The difference between our results and those of Husain-Seahra may have at least two 
origins. One is the diffuculty in dealing numerically with the nonparabolic system 
()3.161 13.17P . Another is the difficulty in numerically modelling fall-off conditions on a 
noncompact manifold using boundary conditions at a fixed boundary. Parabolic equa- 
tions exhibit instantaneous action-at-a-distance, though it is typically exponentially 
suppressed. Our numerical runs appear to produce some positive scalar curvature at 



®The difference appears due to numerical computation of second spatial derivatives, which grow to 
much larger magnitude when tp is the second (nonzero) option in (|3.4[) than when -0 = 0. 
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Figure 5: Area of the minimal surface at r = 1 as a function of flow time for initial data 
give by a time-symmetric slice of an (n + l)-dimensional exterior Tangherlini metric for 
n = 4, 5, 6, and 7. 
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Figure 6: Area of the minimal surface at r = 1 as a function of flow time, for a = 2, 
3, 4, and 3. The DeTurck diffeomorphism is generated by the difference between the 
flowing connection and a flat connection. 
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Figure 7: Scalar curvature near the outer boundary at r = 100 for a = 3, pfotted 
at three different times, corresponding roughly to 10%, 50%, and 90% of the time 
to collapse of the minimal surface. At early times, a disturbance is present near the 
boundary at r = 100, but it appears to dissipate without propagating inward. (For 
comparison, at r = 1, the initial scalar curvature is i?(0, 1) = —4). 

the large r boundary, which must be chosen sufficiently distant to suppress the effect 
on the minimal surface. The effect is shown in figure 7. For the figure, we integrated 
over r £ [1,100] but chose to display only the portion r = [90,100]. At early times 
we see a small pulse of positive scalar curvature concentrated at infinity, presumably 
attributable in some way to the boundary condition. As time progresses, the pulse 
dissipates. The figure shows three times, one shortly after the fiow begins, one at 
around the midpoint in time of the flow, and one just before collapse of the minimal 
surface. We chose to display a = 3 but the behaviour is similar for all a. On the other 
hand, as discussed in subsection 3.2.2, the boundary conditions in [LQl imply a concave 
boundary at r = Tc and, concomitantly, appear to produce a source of negative curva- 
ture there. Then the maximum principle argument of section 2 would not apply. The 
critical behaviour observed in [10] may be a property of boundary conditions, rather 
than evidence for unstable regions in the space of initial data for asymptotically flat 
metrics. 
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